AD-76  9  585 


DESIGN  AND  ANALYSIS  OF  90-DEGREE  PHASE 
DIFFERENCE  NETWORKS 

W  .  W  i  e  b  a  c  h 

Harry  Diamond  Laboratories 
Washington,  'D.  C. 

May  1973 


National  Tsctinical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal- Road,  Springfield  Va.  22151 


\Z».  HtPOrtT  SECURITY  CLASSIFICATION 


UNCLASSIFIED 


Security  O*  ossification  ... 


DOCUMENT  CONTROL  DATA  •  R  S>  D 

(Sacurlty  claaslll*: at  Ion  of  tltU,  body  of  abattact  and  tndaxtng  annotation  mu  at  ba  anlarad  whan  tha  ova/all  rarporf  I*  ela^aKlad 


I  ORIGIN  A  TING  ACTIVITY  (Cotgcraia  author) 

Harry  Diamond  Laboratories 
Washington,  D.  C.  20438 


3.  REPORT  TITLE 

DESIGN  &  ANALYSIS  OF  90  DEGREE  PHASE  DIFFERENCE  NETWORK 


UNCLASSIFIED 


4.  OESCRiPTive  NOTES  (Typ*  of  raport  and  Inctualvo  datoa) 


9.  auTmor<5>  (Hrat  natoa,  tnlddt*  Initial,  laat  ttama) 


W.  Wiebaeh 


*  REPORT  DATE 

June  1973 


ta.  CONTRACT  OR  GRANT  NO- 


76.  NO.  or  REFS 


l»*.  ORIGINATOR*!  REPORT  NUMBER!*) 


b.  pnojtcTMO.  DA-11B15011 


HDL-TM-73-4 


c.  AM  CMS  Code:  4210.18.08102 


46.  OTHER  REPORT  nO(S>  (Any  othat  ntun6«r«  that  may  ba  aaatjnad 
thla  raport) 


«.  HDL  Proj.  No.  47311 


10.  OtSTRI3UT.CN  STATEMENT 

APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


12  SPONSORING  MILITARY  f  CTIVlT** 

U.  S.  Army  Matei'iel  Command 


IS.  ABSTRACT 


The  parameters  of  a  90  degree  phase  difference  (Dome)  filter,  and  their  relationship 
to  the  filter  characteristics,  are  derived  and  discussed.  The  design  of  two  RC  implemen¬ 
tations  and  their  worst-case  analyses  are  shown. 


Peprodoccd  b> 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

U  S  0*pafim»nt  nf  CoRunerc# 
5p'ingf'»M  VA  ??1  V 


473 


REPLACES  OO  FORM  147S.  <  JAN  M,  WHICH  IS 
OBSOLETE  FOR  ARMY  USX. 


UNCLASSIFIED 


UNCLASSIFIED 


"Security  Clectlflcatlcn 


key  no  not 


Dome 

Filter 

Phase 

Difference 

Network 


s 


Security  CUntlflcsUo 


ABSTRACT 


The  parameters  of  a  90~degree  phase  difference  (Dome)  filter,  and  their 
relationship  to  the  filter  characteristics,  a  'e  derived  and  discussed.  The  design 
of  two  RC  implementations  and  their  worst-  ease  analyses  are  shown. 


TABLE  OF  CONTENTS 


Section  Title  Page 

ABSTRACT .  2 

1  i  introduction .  7 

2  OBJECTIVE  OI  THIS  REPORT .  8 

3  DERIVATION  OF  THE  FILTER  PARAMETERS  .  9 

3.1  Frequency  Response  of  che  Phase  Anglo .  9 

3.2  Condition ’or  Equal  Maxima .  12 

3.3  Conditions  for  Equal  Ripple  Around  90  Degrees . .  15 

3.4  Bandwidth .  ]5 

3.5  Filter  Parameters  for  a  Specified  Bandwioth  ........  17 

3.6  Filter  Parameter's  for  a  Specified  Deviation  .  18 

4  RC  REALIZATION  WITH  BALANCED  GENERATOR .  22 

4.1  All-pass  Conditions  . 22 

4.2  Dimensioning . 24 

4.3  Effect  of  Component  Tolerances  .  26 

4.4  Amplifier  Considerations . . . 28 

4.5  Example .  30 

5  RC  REALIZATION  WITH  DIFFERENTIAL  AMPLIFIER .  35 

5.1  All-pass  Condition  . .  39 

5.2  Dimensioning . 37 

5.3  Effect  of  Component  Tolerances  .  38 

5.4  Example .  38 

6  SUMMARY .  30 

LITERATURE  CITED .  43 

APPENDIX  ] .  45 

APPENDIX  2 .  :? 

APPENDIX  3  53 

DISTRIBUTION  LIST  .  55 


Preceding  page  blank 


LIST  OF  ILLUSTRATIONS 


Figure  T'tle  Page 


1  Phase  response  versus  frequency  obtained  with 

I  .mo's  Dimensioning  . . .  7 

2  Equal-ripple  phase  response  . . 7 

3  Poles  and  zeros  of  second-degree  all-pass .  9 

4  Basic  dome  filter  . . 12 

5  Phase  shift  and  pnase  difference  of  an  arbitrary  dome  filter  .  .  13 

6  Frequency  function  z . 13 

7  Definition  of  edge  frequencies  i\t  f 2 .  16 

8  Order  of  calculation,  and  equations  used,  when  the  band 

edges  are  specified . . . .  18 

9  Order  of  calculation,  and  equations  used,  when  the  maximum 

deviation  is  specified  . . . . .  20 

10  Pole  and  zero  location  for  >  2  . .  20 

11  All-pass  circuit  with  balanced  generator  .  21 

12  Possible  circuit  parameters  for  s  =  4 .  25 

13  Dome  filter  fed  by  simple  phase  splitter,  not  recommended  .  .  29 

14  Dome  filter  fed  by  two  transistors .  29 

15  Dome  filter  fed  by  two  operational  amplifiers  ...........  30 

16  Basie  all-yass  network  . 35 

17  Sebol’s  ali-pass  circuit .  35 

18  Solutions  of  r,  k  vs  c  for  s  =  4  37 

LIST  OF  TABLES 

Table  Title  Page 


I  Drive  Circuit  With  Bandwidth  Generator . . .  33 

II  Dome  Effect  With  Differential  Amplifier .  40 

Computer  Program  For  the  Worst  Case  Analysis  of  the 

Dome  Filter  Circuit  With  Balanced  Generator .  46 

Subprogram  AECT  Required  By  Programs  in  Appendixes 

Appendixes  I  and  2  . . . . .  49 

Computer  Program  For  the  Worst  Case  Analysis  of  the 

Dome  Filt  -.r  Circuit  With  Operational  Amplifier .  52 


(3 


1.  INTRODUCTION 


The  "Dome  filter"  was  originally  described  by  R.  B.  Dome1  and  consists  of 
a  pair  of  second-degree  all-pass  networks  with  a  phase  shift  difference  in  the 
vicinity  of  90  degrees  over  a  -vide  bandwidth.  Dome's  article  shows  three  differ¬ 
ent  implementations  -  one  I.C  and  two  RC  circuits  -  and  gives  equations  for  the 
component  values  based  on  an  apparently  arbitrary  (within  limits)  parameter  (s). 

He  suggests  frequency  parameters  for  the  two  all-pass  sections  which  result  in 
a  one-sided  phase  difference  curve  as  shown  in  figure  1. 

Antony  published  a  report2  on  Dome  filter  design  with  a  computer  program 
for  component  dimensionii  g.  He  points  out  that  there  is  a  relationship  between 
the  bandwidth  and  the  ripple  of  the  phase  error,  which  is  determined  by  Dome's 
parameter  s.  He  suggests  to  shift  the  phase  difference  response  down,  as 
shown  in  figure  2,  to  obtain  equal  positive  and  negative  deviations  from  the  nominal 
value.  For  a  given  bandwidth,  this  design  modification  almost  halves  the  error  or, 
conversely  for  a  given  maximum  error,  it  almost  doubles  the  bandwidth. 


A(f> 


Figure  1.  Phase  response  versus 
frequency  obtained  with 
Dome's  dimensioning. 


Figure  2.  Equal-ripple  phase 
response. 


’R.  B.  Dome,  "Wideband  Phase  Shift  Networks,"  Electronics,  Dec  1946,  p  112 
vol.  4. 

^R.  T.  Antony,  "Dome  Filter  Design  and  Analysis  Program, "  HDL  TM-70-26,  1970. 
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Another  RC  realization  of  the  Dome  filter  with  a  differential  amplifier  was 
given  by  Sebol.3  He  use0  as  filter  parameters  tabulated  data  published  by 
Bedrosian,4  which  also  permit  the  design  of  higher  order  phase  difference  net¬ 
works  by  cascading  second-degree  all-pass  sections.  The  phase  response 
corresponds  to  figure  2. 


2.  OBJECTIVE  OF  THIS  REPORT 

Neither  one  of  the  cited  papers  presents  much  of  the  theoretical  background 
of  the  Dome  filter,  the  knowledge  and  understanding  of  which  can  add  signifi¬ 
cantly  to  the  designer's  confidence.  This  is  one  gap  which  this  report  intends 
10  fill.  While  the  higher-order  filters  are  not  accessible  by  elementary  analy¬ 
sis,  the  second  order  (Dome)  filter  is  readily  so. 

The  cited  papers  pi*esenl  various  graphs  showing  the  relationship  between 
the  filter  parameters,  bandwidth  and  ripple;  Antony's  computer  program  calcu¬ 
lates  the  component  values  from  the  phase  deviation.*  This  report  derives  the 
calculation  of  the  filter  parameters  based,  at  the  designer's  choice,  either  on 
the  required  bandwidth  or  on  the  permissible  phase  error.  It  also  shows  the 
flexibility  which  exists  in  dimensioning  component  values. 

Dome  uses  an  arbitrary  capacitance  ratio  in  his  formulas  for  the  compon¬ 
ent  values,  but  deprives  the  designer  of  one  degree  of  freedom.  Antony  uses 
this  capacitance  ratio  in  his  computer  program.  Also,  the  user  of  his  program 
seems  to  have  little  choice  as  to  the  component  determining  the  impedance  level 
of  the  fil.vr.  However,  as  this  paper  shows,  any  component  of  a  filter  section 
may  be  seier  led  for  this  purpose,  and  the  impedance  of  the  two  sections  need 
not  necessarily  be  the  same.  It  is  deemed  important  to  acquaint  the  designer 
with  the  flexibility  avr liable  in  this  circuit  so  that  he  may  realize  it  with  as  low 
an  impedance  spread  as  possible  and  as  many  as  possible  standard  values  for  the 
capacitors,  because  non-standard  capacito.;?  are  expensive  and  hard  to  get. 

The  calculation  of  the  component  values  for  Sebol's  circuit  is  also  summa¬ 
rized  to  make  it  available  to  the  user  in  the  same  frame  of  reference. 

As  this  paper  shows,  the  calculation  method  for  filter  parameters  and  com¬ 
ponents,  even  with  the  described  full  flexibility,  is  not  unduly  complicated, 
requiring  only  a  desk  calculator. 

Finally,  this  report  presents  worst-case  analyses  of  the  effects  of  compon¬ 
ent  mleranees  on  the  phase  difference  for  a  Dome  circuit  and  the  Sebol  circuit. 
Computer  programs  for  this  purpose,  and  representative  sample  outputs  are 
listed. 


3 Sebol,  R. ,  "Design  of  Active  90-Degree  Phase  Difference  Networks," 
HDL-TM-68-18 

4Bedrosian,  S.  D.,  "Normalized  Design  of  90-Degree  Phase  Difference 
Networks,"  IRE  Trans,  on  r,lreuit  Theory,  June  1960. 


♦With  some  error:  in  one  particular  case  the  component  values  from  the 
computer  program  were  up  to  three  percent  off. 
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3.  DERIVATION  OF  THE  FILTER  PARAMETERS 


3.1  Frequency  Response  of  the  Phase  Angle 

The  transfer  function  of  a  second-degree  all-pass  network  with  the  complex 
frequency  p  is  (fig.  3)s 


V,(P)  (P  Pi)  (P  +  Pt) 

Vj  (P  -  Pj)  (P  -  PL) 


Where  p^  =  -<*  +  j/J. 

It  has  two  conjugate-complex  poles  and  two  conjugate-complex  zeros  located 
symmetrically  to  the  poles  with  respect  to  both  the  origin  and  jw-axis.  The 
frequency  response  can  be  written  as 


v2(")  [-  rt  +  j  ("+0)]  [-«  +  j  (“>-/?)] 

Vj  [«  +  j  (»-•/?)]  [«  +  j  )] 

ItL'lll 

~  lV.agn. - -  Magn.  Ix. 

Ifi  * 


(2) 


Figure  3.  Poles  and  zeros  of  second-degree  all— pj  ss. 


see,  for  instance,  L.  Weinberg,  "Network  Analysis  and  Synthesis,"  Me  Craw 
Hill,  19(>2,  p  2*5. 
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with  the  phase  angle  $  -  ^  t  -  <£,.  -  <f>. . 

1  J  i)  4 

The  phase  angles  of  the  individual  terms  are  expressed  by 

tan^.,  -- 

.)  Oi 

tan</>,  -  ~ 

"t  or 

to  4  fj  f  \ 

ton*i  -rrr  tan(5r_  ti) 

tan 4>2  =  — ~Jl  -  tan  fir  -  4>.^j 

thus 


(3) 

(4) 

and  substituting  «  2 7r f 

4  -  aretan 
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Comparing  this  result  with  Dome's  equation 


2s  f  fo  (f2-fo  ) 

'  (f2-fo2)2  -  (s'  fo  f)2’ 


shows  that  Dome's  parameters  s,  f  are  related  tooc ,  P  by 


Using  (G)  to  eliminate  a  and  P  from  (3)  yields 


(5) 


<j>  =  2  arctan  —•  •  (8) 


The  function  v  is  often  used  to  describe  the  response  of  tuned  cinuiti?  or 
band  pass  filters  which  are  symmetrical  (although  not  linearly)  abotu  a  fre¬ 
quency  fQ,  corresponding  to  v  =•  0.  So  the  phase  response  of  this  al'-psse  net¬ 
work  is  in  the  same  way  symmetrical  about  f  ,  which  could  be  called  ;i* 
center  frequency;^  (f0)  -  -180°.  0 
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The  dome  filter  usee  two  all-pass  sections  connected  as  in  figure  4,  Of 
interest  is  the  phase  difference  of  the  output  signals: 


- 

2  ' 


arctan 


1  + 


S1S2 

vlv2 


(9) 


where  and  are  defined  by  (7)  with  the  center  frequencies  fol  and  f02» 
respectively.  An  example  of  the  phase  shift  with  some  arbitrary  parameters 
is  shown  in  figure  5. 

3.2  Condition  For  Equal  Maxima 

Setting  the  deri  native  equal  to  zero,  solving  for  the  locations  of  the 

extrema  fmaxl>  ^max2*  *min  and  lettinc.  =  ^max2  would  result  in  a 

condition  between  the  filter  parameters  tor  equal  maxima.  However,  this  is 
too  unwieldy  to  be  done  in  a  general  way.  But  it  turns  out  that 


s 


1 


=  s 


(10) 


is  one  condition  resulting  in  equal  maxima,  as  will  be  shewn.  For  this  purpose 
the  following  new  frequency  function  is  introduced  (fig.  6): 


z 


f  4fm 
fm  f 


(ID 


V  z*1 

vo! 


Figure  4.  Base  dome  filter. 
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fm  fm  fm  2  fm  4  fm  logf 
4  2 


Figure  5.  Phase  shift  and  phase  Figure  6.  Frequency  function  z. 

difference  of  an  arbitrary 
dome  filter 

(foi  =  °-6;  si  =  3-5; 
fo2  =  2'°;  s2  "  4-6> 

where 


tm  =  ■ 

With 

b  = 


(12) 


(13) 


the  center  frequencies  of  the  two  all-pass  sections  can  be  expressed  as 


l 


u 


fm 

b 


b  fm  . 


(14) 
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Using  these  frequency  functions  and  equation  10,  equation  9  can  be  written 
as  a  function  of  z: 


Under  the  condition  of  (10),  the  phase  difference  &<f>  is  obviously  symmetrical 
about  fm,  since  there  are  always  two  values  of  f  (geometrically  symmetrical 
to  frn)  which  give  the  same  z,  and  is  a  function  of  z  only.  So  the  z-vaiue 
Zmax  where  (15)  has  a  maximum  translates  into  two  frequencies,  located  on 
opposite  sides  of  fm,  where &<f>  has  the  same  maximum  amplitude  A^(zmax). 

J  4  I 

zmax  cou^cl  be  f°und  by  setting  =  0,  but  this  procedure  does  not  also 

yield  the  minimum.  The  locaticn  of  the  minimum  can  be  ascertained  by 
re-introducing  the  frequency  into  (15)  and  differentiating  with  respect  to  f.  As 
it  turns  out, 


f  . 
mm 


=••  f 

m 


corresponding  to  zmjn  =  2;  at  this  point,  a^( z)  has  an  absolute,  not  a  relative 
minimum.  The  frequencies  of  the  maxima  are 


f 

max 

f 

m 


^HN±mE 


(16) 


The  values  of  the  function  at  the  minimum  and  at  the  maxima  are  as  follows: 


A.f nriin  _  arclan 


^max  =  arctan 
2 


(17) 


(18) 
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3.3  Conditions  for  Equal  Ripple  Around  90  Degrees 

For  most  practical  applications,  a  phase  difference  of  90  degrees  is 
required.  The  condition  for  equal  positive  and  negative  maximum  errors  with 
respect  to  90  degrees  (fig.  7)  is 


f>max  -♦  A<j 
2 


■:  90° 


\  2  2  /  i-an-a-^S5lan^-‘i 


-  =  tan  90  -oo 


1  -  (an  i|2S2  tt„  ifUS.  =  0 


which,  with  (17)  and  (18),  leads  to 


This  equation  may  be  combined  with  one  other  condition  to  solve  for  the 
parameters  s  and  b:  either  the  maximum  deviation  or  the  required  bandwidth. 

3.4  Bandwidth 

The  edge  frequencies  (fig.  7)  that  describe  the  useful  frequency  range  are 
defined  as  the  frequencies  where  the  error  of  the  phase  difference  is  the  same 
as  at  the  minimum,  6,  which  is  equal  in  magnitude  to  the  deviation  at  the 
maximum.  At  this  band  edge,  the  frequency  function  z  has  the  value 


il_  ,  fm  ^2_  ,  fm 
ZL  fm  fj  =  fm  T  f2  5 
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solved  for  f  and  zT  : 

m  ij 


z 


L 


or  solved  for  fj  and  f0 


Figure  7.  Definition  of  edge  frequencies  f 
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It  may  be  solved  numerically  using  Newton’s  method: 


s2  "  sl' 


the  derivative  is 


r(si) 

F'(si) 


F’(s)  ■*  8D6s?  -  96D2s5  +  128s3  (  D4  +2D2-^  +  512s. 


After  the  solution  of  s  has  been  found,  b  is  calculated  from  (27).  The 
resulting  .bandwidth  is  obtained  again  by  considering  that,  at  the  band  edge 


z  =  zL, 


A^(Zj  ^  =  4<£min  . 


Substituting  both  parts  of  this  equation  by  (15)  and  (17)  and  solving  for  Zjj 


7. 


L 


•iiy-M)*] 


(29) 


The  edge  frequencies  are  obtained  from  (22).  Figure  9  summarizes  the  order 
of  calculations.  The  following  tabulation  shows  some  numerical  examples: 


s 

2  deg 

3  deg 

4  deg 

!  8 

3.S68144 

4. 102830 

4.325051 

b 

2.050514 

2.110651 

r,.  165476 

ZL 

4.260048 

5.076946 

5.901764 

1 

*— 1 
<~4 

1 

_ 

16.0358 

. 

23.7332 

32.8003 

if! 

) 

% 
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Figure  9.  Order  of  calculation,  and  equations  used, 
when  the  maximum  deviation  is  specified. 


For  8  •-  0,  the  smallest  possible  value  of  s  is  obtained  (s  =  3.107543).  Going 
back  to  the  original  parameters  « ,  p  of  the  all-pass  network  and  solving  (6) 
tor  /?: 


shows  that,  for  all  practical  filters,  p  is  imaginary,  j P  is  real  and  ali  poles 
and  zeros  are  located  on  the  real  axis  (fig.  10).  Hence,  the  network  s.  suitec 
for  RC  realization. 


Figure  10.  Pole  and  zero  location  for  s  >  2. 


4.  RC  REALIZATION  WITH  BALANCED  GENERATOR 

Rr  n«fe  of  ^  Impiementeiions  of  the  all-pass  filter  suggested  by  Dome  is  an 

^r?'!?’k  fed  by-  tw°  out'of"Piias£i  signals,  figure  11.  in  view  of  the  -vorst- 
case  analysis  performed  later,  a  possible  difference  of  the  maani!  des'nf  tho 
S'f  *“•  b-n  considered.  The  frequency  relpo 


^  -  (l  +  -iV)C, 

'1  =  c2  +  S 


~  I 

Ocu)2  -  jo,  r — i — _ i _i  i  i 

Rici  “  vd +  vw; 

JL  +  J_  n  i  l” 

M)2  +  j.  iUi ,  _jl.  +  i  _SiE^ 

j_C2  +  c:>  B1(VC3)j+R^c^ 


1  A 11- pass  Conditions 

'Pt,  aH_Pf,SS  conditions  of  circuit  of  figure  11  are  derived  for  \V  =  n 
rhe  general  all-pass  transfer  function  (2)  may  be  written  as 


2  (jw)2  -  2«  j<o  +  <*“  +/?“ 

V  ^  ^ '  r  *  1  p 1  A" 

1  (juj)  +  2«j<4)  +  or  + 


Compared  to  (30),  it  follows  that  (31)  describes  an  all-pass  if 


k**»  ,  i 


C2  +  C3  R1C1  W3 


R1C2  '  »XCl  '  r2C2 


fr 

*.* 


4.2  Dimensioning 

As  shown  by  (34)  the  ti’ansfer  function  of  the  all-pass  circuit  is  expressed  by 
the  angular  center  frequency  a>0  and  the  two  circuit  parameters  c  and  r.  The 
desired  filter  parameter  s  (par.  3.5  and  3.6)  provides  one  condition  between 
then,  so  one  parameter  is  selectable.  Recalling  from  (6)  that 


and  comparing  (34)  with  (31)  shows  that 


it  is  convenient  to  select  the  capacitance  ratio,  this  equation  is  solved  for 
the  resistance  ratio 


VF  = 


-s 

2  vc 


(35) 


The  dimensioning  continues  by  calculating  the  ratios  R3/R2  and  C2/C3,  which 
are  equal  according  to  (32)  and  follow  from  c,  r  using  (33); 


The  ratio  c  has  to  be  selected  so  that  both  r  and  a  are  positive  real  numbers 
Figure  12  shows  as  an  example  the  range  of  possible  values  for  the  filter 
parameter  s  =  4. 


Dome  and  Antony  arbitrarily  set 


'  ix&yVFf&Sf, 


which  means 


and 

c  =  s  +  2 


After  a  set  of  circuit  parameters  c,  r,  a  has  been  found,  continue  by 
selecting  one  of  the  three 


resistors  or  capacitors 


With 

nominal  parameters  uq,  a,  c,  .r; 

nominal  components  R^,  Ii2o,  U3o,  C1q,  C2o,  C.jo,  aV  =  0; 
actual  components  P, ,  R„,  R  ,  C,,  C’  CL,  aV; 

i.  It  t)  X  £•  «> 


2C> 


the  transfer  function  of  one  all-pass  section  (31)  can  be  written  as: 


v/cv 


a  R, 


R 


1C 

R1 


C10  1  fc  _ R1 _ 

Cj  ^/c?  \r  C2  C3 

C„„  4  a  C'„ 
20  30 


Equation  (38)  contains  only  the  nominal  circuit  parameters  and  ratios  of  nomi¬ 
nal  to  actual  component  values.  In  this  way,  the  effect  of  percent  ile  variations 
of  components  on  the  frequency  response  can  be  studied  for  a  .more  general 
network.  A  computer  program,  listed  in  appendix  1,  has  been  written  to  caleu- 
le  e  the  worst-ease  phase  difference  for  combination  of  specified  component 
tolerances. 

Tolerances  may  be  specified  for  each  of  the  three  resistors,  for  each  of 
the  three  capacitors  and  for  the  symmetry  of  the  input  voltages.  The  program 
assumes  the  same  tolerances  for  corresponding  elements  in  the  two  channels. 

It  further  uses  the  fact  that  the  toiorance  of  a  pair  of  components  (e.g. ,  R2  in 
channel  1  and  2)  affects  the  response  most  if  the  two  components  are  at 
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opposite  ends  of  their  tolerance  range  (channel  Is  R2  =  R20  +  A^2* 
channel  2:  R2  r-  R20  “  AR2).  The  program  goes  through  all  possible  combina¬ 
tions  of  tolerances  {21  combinations  if  all  tolerances  are.  specified  as  non-zero); 
and  in  each  case  sweeps  through  the  frequency  range  and  notes  the  phase  differ¬ 
ences  at  the  extrema  and  the  band  edges.  The  smallest  and  largest  values  for 
all  combinations  at  each  of  these  points  is  printed  out. 

The  band  edges  are  specified  as  the  ratio  fg/fj,  from  which  the,  filtei 
parameters  s,b,  s  are  immediately  computed  (eq.  23,  24,  25).  The  running 

f 

frequency  variable  is  x  ■=  —  ,  and  it  is  swept  from 

rn 


The  frequency  variable  actually  substituted  into  equation  (38)  is  (using  (12) 
and  (13)): 


For  Channel  1 


For  Channel  2 


~  X  .  I  . 
">  2  ~  fm  b  ’ 


Since  the  absolute  frequency  limits  need  not  be  specified,  the  worst-case 
analysis  is  valid  for  any  filter  of  such  configuration  characterized  by  band  width 
f2/i:  and  implementation  parameters  c(chan.  1),  c(chan.  2). 

4.4  Amplifier  Considerations 

The  Dome  filter  implementation  treated  in  thi  *  chapter  must  be  fedbv  t»,vo 
signals  of  equal  amplitude  and  opposite  phase  (fig.  11).  A  simple  one-transi.  '*or 
phase  splitter  as  in  figure  13  probably  does  not  suffice  when  a  phase  error  of  a 
few  degrees  is  desired,  *  The  change  of  the  emitter  load  C2H  R2  with  frequency 
reflects  onto  the  equivalent  collector  voltage  source. 

This  interaction  is  considerably  reduced  in  the  t.wo-transistor  circuit  of 
figure  14.  It  should  work  satisfactorily  when  the  emitter  load  impedances  Ri 
and  Re  II  R2  ||  C2,  multiplied  by  the  transistor  current  gain,  each  are  large 
compared  to  the  generator  impedance  Z^. 


♦supported,  by  experience  relayed  by  K.  Sann,  HDL 


Impedance  relationships  need  not  be  considered  when  feeding  the  Dome  filter 
by  a  pair  of  intcgrated-cireuit  operational  amplifiers  (fig.  15).  One  operational 
amplifier  works  in  the  inverting  mode,  the  other  in  the  noninverting  mode. 

With  unity  gain,  as  shown,  the  voltage  symmetry  depends  upon  the  two  resis- 

a  R  f 

tors  Rp  being  equal.  With  a  resistor  tolerance  —  ,  the  gain  tolerance  is 


AV  =  2 


The  amplifiers  may  provide  voltage  gain;  in  that  case  there  are  two  pairs  of 
resistors  whose  ratio  must  be  equal. 

A  phase  difference  between  the  outputs  of  the  inverting  and  noninverting 
amplifier  is  possible,  particularly  at  frequencies  approaching  the  upper  limit  of 
the  operational  amplifier;  this  error  is  not  included  in  the  worst-case  analysis. 

4.5  Example 

Design  a  phase  difference  network  of  90  degrees  with  a  maximum  nominal 
deviation  of  ±2  degrees  and  a  bandwidth  as  great  as  possible,  extending,  if 
pos'  ible,  from  600  to  10000  Hz. 

8-2° 


Figure  15.  Dome  filter  fed  by  two  operational  amplifiers. 


*»ra»PW'«^ 


According  to  the  scheme  of  figure  9: 

s  ■=  3.868144 
b  =  2.050534 
Zj  «  4.260048 

Using  the  desired  edge  frequencies  to  find  a  suitable  center  frequency 

f  « VGOO'IOOOO  Hz  =  2449.5  Hz - ► 

m 

f  «-  2450  Hz. 
m 

The  actual  band  limits  are  then 

fx  =  610.9  Hz 
f2  =  9826  Hz 

and  the  other  characteristic  frequencies  are 

f  .  *  1194.82  Hz 
ol 

f  .  -  5023.76  Hz 
o2 

f  ,  *  971. 3  Hz 
maxi 

f 

ma.\2 


6180  Hz 


It  is  generally  advantageous  to  select  capacitor,  rather  thaji  resistor 
values,  to  be  standard  values.*  For  each  all-pass  section,  two  capacitors 
can  be  selected.  In  this  example,  and  C2  were  chosen  by  trial  and  error  so 
that  C3  comes  close  to  a  standard  value  too.  Depending  on  the  permissible 
actual  phase  error,  this  standard  value  may  be  usable.  If  nothing  else,  it 
makes  breadboarding  easier. 


Section  1 

Section  2 

C  12  nF 

Cx  =  6.8  nF 

C2  ---  2.  2  nF 

C2  =  1  nF 

c  =  5.454545 

c  =  6.800000 

r  =  0.158105 

r  a  0.196085 

a  =  0.464393 

a  =  0.457115 

C3  =•  4.737  nF 

C3  -  2.138  nF 

R  -  10.31  kn 

Rx  =  5.380  kn 

R0  =  65.20  k« 

R2  *  27.44  kn 

R3  ^  30.28  kO 

R3  =  12.54  kn 

The  results  of  the  worst-case  analysis  are  shown  in  table  I.  After  listing 
of  the  parameters,  the  frequency  is  swept  once  with  the  nominal  values  to  ensure 
that  the  input  data  are  correct.  The  next  seven  results  show  the  effect  of  the 
variation  of  one  component  at  a  time.  As  it  is,  Rj  and  C2  affect  the  phase  dif¬ 
ference  mainly  at  the  higher  frequencies,  II2  and  mainly  at  the  lower  frequen¬ 
cies.  A  voltage  unbalance  shows  up,  although  to  a  lesser  degree,  at  both  ends 
of  the  band,  'flic  effect  of  variations  of  R3  and  C3  is  only  about  one-fourth  of  that 
of  the  other  resistors  and  capacitors. 

With  one  percent  tolerance  of  ail  components  (except  one-half  percent  for 
Rjr  if  the  circuit  shown  in  figure  15  is  used),  the  possible  phase  error  is  almost 
three  degree,  in  addition  to  the  nominal  two  degree  error  with  ideal  compon¬ 
ents.  If  all  resistor  tolerances  are  reduced  to  one-fourth  percent,  the  possible 
error  is  close  to  four  degrees,  or  twice  the  nominal  error.  If  standard  values 
(±1  percent)  are  ordered  for  Cj,  its  possible  variation  may  be  up  to  1.7  percent 
from  the  theoretical  value;  as  the  last  output  shows,  the  error  is  increased  only 
very  slightly. 

‘Although  one-percent  capacitors  may  be  ordered  with  any  arbitrary  nominal  value, 
nonstandard  values  have  higher  prices  and  longer  delivery  time.  This  was  found 
(o  be  the  case  with  ceramic  capacitors  from  Aerovox,  Eric  and  Vitramon. 
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TABLE  I.  DOME  CIRCUIT  WITH  BALANCED  GENERATOR 


NOMINAL  ?*ND«I0TH  F 2/FI*  16.0800 

NOM.  DERIVED  PARAM.  S*3. 867941  8*2.050461  DELTA*2.00 


MOM.  REAUZ.  CHAN.  1  A*  0.46434  C*  5.45450  R  =  0.15811 

CHAN. 2  0.45706  6.80000  0.19610 

*********»*********'****>************************♦********* 

FI  1 . MAX  MIN  2. MAX  F2 


OP  NOM. /DEG 


88.00  92.00 


COMPONENT  VARIATION  U/-PCT) 
(OPPOSITE  IN  CHAN. 1.2) 


OP  MAX  /OEG 
DP  MIN  /OEG 


88.14  92.26 

87.86  91.74 


COMPONENT  VARIATION  (+/-PCT) 
(OPPOSITE  IN  CHAN. 1,2) 


OP  MAX  /OEG 
OP  MIN  /OEG 


88.71  92.68 

87.29  91.32 


COMPONENT  VARIATION  (♦/-PCT) 
(OPPOSITE  IN  CHAN. 1,2) 


OP  MAX  /OEG 
OP  MIN  /OEG 


83.25  92.25 

87.75  91.74 


COMPONENT  VARIATION  (+/-PCT) 
(OPPOSITE  IN  CHAN. 1,2) 


OP  MAX  /OEG 
DP  MIN  /DEG 


88.96  92.93 

87.04  91.06 


COMPONENT  VARIATION  (♦/-PCT) 
(OPPOSITE  IN  CHAN. 1,2) 


DP  MAX  /CEG 
DP  MIN  /OEG 


88.02  92.12 

87.98  91.88 


COMPONENT  VARIATION  (♦/-PCT) 
(OPPOSITE  IN  CHAN. 1,2) 


DP  MAX  /OEG 
OP  MIN  /OEG 


88.12  92.14 

87.88  91.85 


COMPONENT  VARIATION  (♦/-PCT) 
(OPPOSITE  IN  CHAN. 1,2) 


83.62  92.51 

97.37  91.48 


MI 

N 

2. MAX 

F2 

XSSX2 

*  :sis 

EXXKSSX 

88. 

00 

92.00 

88.00 

Rl 

R2 

R3 

C 1 

C2 

C3 

DV 

1.00 

-0.00 

-0.00 

-0.00  -0.00 

-0.00 

-0.00 

88. 

64 

92.93 

86.9  b 

87. 

35 

91.06 

87.04 

Rl 

R2 

R3 

C) 

C2 

C3 

OV 

0.00 

1.00 

-0.00 

-0.00  -0.00 

-0.00 

-0.00 

88. 

43 

92.11 

88. 'V 

87, 

57 

91.89 

87.  78 

Rl 

R2 

R3 

Cl 

C2 

C3 

OV 

■O.OG 

-0.00 

1.00 

-0.00  -0.00 

-0.00 

-o.co 

88. 

C  1 

92.15 

88.12 

87. 

79 

91.85 

87.68 

Rl 

R2 

R3 

Cl 

C2 

C3 

OV 

•0.00 

-C.00 

-0.00 

1.00  -0.00 

-0.00 

-0.00 

88. 

65 

92.26 

86  » 14 

87. 

35 

91.74 

87.86 

Rl 

R2 

R3 

Cl 

C2 

C3 

OV 

o 

o 

• 

o 

-0.00 

-0.00 

-0.00 

1.00 

-0.00 

-o.co 

68. 

43 

92.67 

88.70 

87. 

57 

91.32 

87.30 

Rl 

R2 

R3 

Cl 

C2 

C3 

OV 

■0.00 

-0.00 

i  -0.00 

-0.00  -0.00 

1.00 

-0.00 

88. 

21 

92.26 

88.26 

87. 

79 

91.74 

87.74 

Rl 

R2 

K3 

Cl 

C2 

C3 

OV 

'0.00 

-0.00  -0.00 

-0.00  • 

-0.00 

-0.00 

1.00 

87. 

99 

92.51 

*8.62 

87. 

99 

91.48 

87.37 

OP  MAX  /OEG 
OP  MIN  /OEG 


TABLE  I.  DOME  CIRCUIT  WITH  BALANCED  GENERATOR  (Continued) 


COMPONENT  VARIATION  U/-PCT)  R1  R2  R3  Cl  C2  C3  OV 

(OPPOSITE  IN  Cf'AN.  1,2-  - 

1.00  1.00  1.00  1.00  1.00  1.00  1.00 

OP  MAX  /DEG  90.80  94.89  90.58  94.89  90.80 

QP  MIN  /DEG  85.15  89.07  85.41  89.07  85.15 

COMPONENT  VARIATION  (+/-PCT)  R1  R2  R3  Cl  C2  C3  OV 

(OPPOSITE  IN  CHAN. 1,2)  - - - 

0.25  0.25  1.00  1.00  1.00  1.00  0.50 

OP  MAX  /OEG  89.87  93.94  89.78  93.86  89.77 

OP  MIN  /OEG  86.11  90.05  86.22  9G.13  86.21 

COMPONENT  VARIATION  (♦/-PCT)  R1  R2  R3  Cl  C2  C3  OV 

(OPPOSITE  IN  CHAN. I, 2)  - 

0.25  0.?c  1.00  1.00  1.00  1.70  0,50 

OP  MAX  /OEG  89.95  94.04  89. 9.  94.04  89.95 

OP  UN  /OEG  86.03  89.95  86.07  C9.95  86.03 
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5.  RC  REALIZATION  WITH  DIFFERENTIAL  AMPLIFIER 

While  the  circuit  discussed  in  the  previous  section  may  be  fed  by  operational 
amplifiers,  the  implementation  suggested  by  SEI30L  needs  an  operational  ampli¬ 
fier  as  an  active  part.  It  is  based  on  the  lattic  network  of  figure  16  with  one 
terminal  of  the  input  voltage  grounded  and  the  output  floating.  In  the  actual  cir¬ 
cuit  of  figure  17,  this  output  is  connected  to  a  differential  (operational)  amplifier, 
and  the  resistive  divider  is  also  utilized  as  the  feedback  path.  The  transfer  func¬ 
tion  is  slightly  different  from  the  one  of  figure  16  and  the  voltage  gain  is  unity. 

The  frequency  response  of  the  circuit  of  figure  17  is 


Figure  16.  Basic  all-pass  Figure  17.  oebol's  all-pass  circuit, 

network. 
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5.1  All-pass  Condition 

Comparing  equation  (39)  with  the  general  all-pass  transfer  function  (31)  shows 
that  the  constant  terms  of  numerator  and  denominator  are  always  identical,  and, 
comparing  with  (16) 


o  **■  -  ■  ~~  * 

°  Jc1c9r1r2 


'file  other  condition  for  an  all-pass  function  is  then 


4  '  C2R2  +  'ClkR2  =  ClRi  +  SR2  +  C1R2 


or  with  the  abbreviations 


13 


the  all-pass  condition  is 


k  -  -^2ci  1. 
r 


If  it  is  met,  the  transfer  function  may  be  written  as 


9 


'""■*75SW?Wi jLr’^W&fi 


This  all-pass  realization  contains  only  two  capacitors,  both  of  which  can  be 
selected.  Furthermore,  they  may  be  the  same.  These  are  two  attractive  features 
of  this  circuit. 

After  C^,  C2  and  thus  c  have  been  selected,  r  and  k  are  found  from  equa¬ 
tions  (42)  and  (40); 


R 


2 


U> 

O 


VciV 


* 


and  follows  from  r.  The  absolute  values  of  R3,  R4  do  not  enter;  their  ratio 
is  given  by  k.  The  impedance  level  should  be  between  the  operational  amplifier 
input  and  output  impedances. 

5.3  Effect  of  Component  Tolerances 

A  worst-case  analysis  has  been  worked  out  similar  to  the  one  described  in 
section  4.3.  The  compuier  program  is  listed  in  appendix  2.  The  analysis 
considers  variations  of  the  passive  components  only,  not  imperfections  of  the 
operational  amplifier.  As  the  latter  is  involved  very  tightly  in  the  performance 
of  the  circuit,  and  finite  gain,  input  impedance  or  common  mode  rejection  cer¬ 
tainly  show  an  effect,  the  analysis  is  not  as  useful  as  in  the  other  case.  It  shows 
necessary  tolerances  for  the  passive  components;  however,  keeping  these  toler¬ 
ances  does  not  in  itself  guarantee  the  desired  performance. 


5.4  Example 

Requirements  and  filter  parameters  as  in  section  4.5.  For  each  channel, 
both  capacitors  are  selected,  and  selected  equal.  For  r,  the  +  sign  in  equa¬ 
tion  (42)  is  used,  because  it  results  in  a  lower  impedance  spread. 


Section  1 

Section  2 

c  -  1 

c  =  1 

*  C0  -  10  nF 

Cj  -  ,C2  --  2. 2  nF 

r  =  2.646158 

r  =  2.646158 

k  3.755813 

k  =  3.7558J3 

Rg  8. 19  kfi 

It  *•  8.  85  kfi 

ilj  21.7  kfi 

Rj  =  23.4  kfi 

R.j  5.6  kfi 

R.,  5. 6  kfi 

ii  j  '  21.0  kfi 

R.  «  21.0  kfi 

4 
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The  result  of  the  worst-case  analysis  is  shown  in  table  II.  R^  and  C2 
affect  the  phase  difference  mostly  at  lower  frequencies,  1*2  and  mostly  at 
higher  frequencies.  Deviations  of  R;j  and  R4  have  a  lesser  effect  at  both  ends 
of  the  band,  none  in  the  center.  A  tolerance  of  one  percent  of  all  passive  com¬ 
ponents  results  in  a  worst-case  error  of  almost  3.5  degrees  in  addition  to  the 
nominal  error  of  two  degrees.  Resistor  tolerances  of  0.  25  percent  approxi¬ 
mately  halve  the  excess  error. 

The  last  result  of  table  II  shows  that  somewhat  different  circuit  parameters 
give  about  the  same  tolerance  sensitivity. 

The  sensitivity  of  this  implementation  is  also  about  the  same  as  the  one 
treated  in  section  4. 


6.  SUMMARY 

This  report  describes  a  simple  method  of  determining  the  parameter  and 
component  values  and  tolerances  for  a  class  of  90  degree  phase  difference 
networks.  After  a  short  discussion  of  previously  published  papers  on  this  sub¬ 
ject  and  the  limitations  of  their  theoretical  treatment,  the  transfer  function  on 
second  order  all-pass  networks,  the  parameters  of  two  such  networks  providing 
a  90-degree  phase  difference  between  their  output  signals,  and  the  relationships 
of  these  parameters  to  the  phase  ripple  and  bandwidth  are  derived.  Simple  and 
straightforward  computational  methods,  requiring  only  a  desk  calculator,  for 
determining  the  filter  parameters  based  on  either  the  required  bandwidth  or  the 
permissible  phase  deviation  are  then  outlined. 

Two  different  inductorless  networks  for  the  implementation  of  the  all-pass 
function  and  their  respective  advantages  and  disadvantages  are  described,  and 
equations  defining  component  values  are  derived.  Here,  as  with  the  filter 
parameters  above,  the  flexibility  and  options  available  to  the  designer  are 
stressed.  For  each  circuit,  a  computer  program  is  listed  and  discussed  which 
permits  the  calculation  of  the  worst  effect  of  component  tolerances  on  the  phase 
difference.  Transistor  and  integrated-circuit  amplifiers  suitable  for  the  active 
parts  of  the  circuits  are  considered  briefly.  Computational  methods  are  illus¬ 
trated  by  examples  where  the  filter  parameters,  the  circuit  components,  and 
worst-case  phase  errors  for  both  implementations  are  calculated. 
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TABLE  II.  DOME  CIRCUIT  WITH  DIFFERENTIAL  AMPLIFIER 

2.  DEG.  EQUAL-RIPPLE  90  DEG.  DIFF.  NW.  WORST  EFFECT 

OF  COMP.  VARIATIONS 


NOMINAL  BANDWIDTH  F2/F1»  16.0800 
NON.  DERIVED  PARAM.  S*3. 867941 


8*2*050461  0ELTA*2.00 


N0M.  REAlIZ.  CHAN. I  C*  1. 00000  K*  3.75593  R=  2.64575 

CHAN. 2  1.00000  3.75593  2.64575 

**********************  I********************-***  ****  ******** 

FI  l.NAX  MIN  2. MAX  F2 

OP  NOM./OEG  88.00  92.00 

CGMPONENT  VARIATION  (+/-PCT)  R1 
(Ot / DSI TE  IN  CHAN. 1.2) 


DP  MAX  /DEG 
DP  MIN  /DEG 


88.96 

87.04 


92.93 

91.06 


COMPONENT  VARIATION  (♦/-PCT)  R1 
(OPPOSITE  IN  CHAN. 1.2) 


DP  MAX  /DEG 
DP  H!N  /DEG 


88.14 

87.86 


92.26 

91.74 


COMPONENT  VARIATION  (+/-PCT)  R1 
(OPPOSITE  IN  CHAN. 1,2) 


OP  MAX  /OEG 
DP  MIN  /DEG 


88.14 

87.86 


92.26 

91.74 


COMPONENT  VARIATION  (♦/-PCT)  R1 
(OPPOSITE  IN  CHAN. 1,2) 


DP  MAX  /OEG 
DP  MIN  /DEG 


88.96 

87.04 


92.93 

91.06 


COMPONENT  VARIATION  (+/-PCT)  Ri 
(OPPOSITE  IN  CHAN. I, 2) 


DP  MAX  /DEG 
DP  MIN  /OEG 


88.65 

87.35 


92.53 

91,46 


COMPONENT  VARIATION  (♦/-PCT)  Rl 
(OPPOSITE  IN  CHAN  1,2) 


DP  MAX  /DEG 
OP  MIN  /DEG 


88.65 

87.35 


92.54 

91.47 


COMPONENT  VARIATION  (♦/-PCT)  Rl 
(OPPOSITE  IN  CHAN. 1,2) 


OP  MAX  /OEG 
OP  HIM  /OEG 


91.45 

84.45 


95.47 

88.51 


88. 

00 

92.00 

88.00 

Rl 

R2 

Cl 

C2 

R3 

R4 

1.00 

-0.00 

-0.00 

-0.00  - 

0.00 

-0.00 

88. 

65 

92.26 

86.14 

87. 

35 

91.74 

87.86 

Rl 

R2 

Cl 

C2 

R3 

R4 

-0.00 

1.00 

-0.00 

-0.00  - 

0.00 

-0.00 

68. 

64 

92.93 

88.96 

87. 

35 

91.06 

87.04 

Rl 

R2 

Cl 

C2 

R3 

R4 

-0.00 

-c.oo 

1.00 

-0.00  -0.00 

-0.00 

88. 

64 

92.93 

88.96 

87. 

35 

91.06 

87.04 

Rl 

(T' 

Cl 

C2 

R3 

R4 

-0.00 

-0.00 

i  -0.00 

1.00  - 

■0.00 

-0.00 

88. 

65 

92.26 

88.14 

87. 

35 

91.74 

87.86 

Rl 

R2 

Cl 

C2 

R3 

R4 

-0.00 

-0.00 

i  -0.00 

-0.00 

1.00 

-0.00 

87. 

99 

92.53 

88.65 

i  87. 

99 

91.46 

87.35 

Rl 

R2 

Cl 

C2 

R3 

R4 

-0.00 

-0.00  -0.00 

-0.00  -0.00 

1.00 

88. 

00 

92.54 

88.65 

'  88. 

00 

91.47 

87.35 

Rl 

R2 

Cl 

C2 

R3 

R4 

1.00 

1.00  1.00 

1.00 

1.00 

1.00 

1  90. 

59 

95.43 

91.45 

85. 

39 

88.51 

84.46 

*&ie{*ffVfg&>&f>'&VK  (wa»‘  m^v.s#  Wi  ^53T’j!?ffr!BP5^se  qjeisFg^*!/i!!Z^TP<>i3q/w&*tmj!a}%8mmK*Fftr®^mi&*fiSvrqgsw 


fc 

f 

if 

< 


table  11.  DOME  CIRCUIT  WITH  DIFFERENTIAL  AMPLIFIER  (Continued) 

COMPONENT  VARIATION  I+/-PC 
(OPPOSITE  IN  CHAN. 1,2) 


OP  MAX  /0E(r 
OP  MIN  /DEG 


89.69 

86. 2< 


NON.  REAL 1Z.  CHAN.l 
CHAN. 2 

FI 

OP  NOM./CEG  88.00** 

COMPONENT  VARIATION  (♦ 
(OPPOSITE  IN  CHAN. 1,2) 


93. T6 
90.23 


C*  0.50000 
0.50000 


Rl  R2  Cl 

C2  R3 

R4 

0.25  0.25  1.00 

1.00  0.25 

0.25 

89.62  93.76 

39.69 

86.39  90.23 

86.29 

K*  3.15261 
3.15261 


R*  1.73519 
1.73519 


l.MAX 

***»«», 

92.00 


OP  MAX  /OEG 
OF  MIN  /OEG 


91.40 

84.51 


95.39 

88.55 


MIN 

2. MAX. 

F2 

Ktimim,,,!,, 

r**»  =  aa 

88.00 

92.00 

88.00 

Rl  R2 

Cl 

C2  R3  R4 

.00  ).00  i.OO 

1.00  1.00  1.00 

90.59 

95.39 

91.41 

85.39 

88.55 

84.51 

•II 


**■** 
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APPENDIX  i 


COMPUTER  PROGRAM  FOR  THE  WORST  -CASE  ANALYSIS  OF  THE  DOME  FILTER 
CIRCUIT  WITH  BALANCED  GENERATOR. 

INPUT  DATA  : 


"As  many  sets  of  data  cards  may  be  used  as  desii’ed. 
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-::r^jj>rK^Jnn^:  ’:  ■ 


-jsx&szi-z 


-w^MBPWSflIW!5Bra*s^?i»» ' 


PROGRAM  LISTING 


OIPENSIQN  A(21,C(2>,R<21* lAm,DC  81  *  Y(  7,21  ,E<4,21  ,PEP(51  ,PEM(51 , 
IP£( 51 

EQUIVALENCE  (  IA(ll,  1 11,  (I  A<  2 ) , 12)  ,  ( I  A(  31 » 13) ,  (IA(4) ,  14) , 
l(IA(51«I51*(IA(61,I61,(IA(7},I71 
PI*3. 1415926 
READ(5, 1)  F 

1  FORMAT! 2F 10. 6) 

Xl*SORT(l./Fl 

X2*$QRT(F} 

G-( (Xi+X<2l**2-4. 1/(X1+X2-2.1 
S«SQRT(G+SvRT1G*G+2.*G*$QkT<  2.*G-4.  1 ) 1 
B-SQRT ( l.+S*S/2.~G+5QRT ( ( l.+S*S/2.-G)**2-i» 1 1 

DELTA«2.*ATAN(<S*S-2.*S4<B-l./B)-(B-l./B  1**21 /(S4S+2MSMB-1./81 
l-(8-l./81**2) 1*180. /PI 
WRITE (6,2  IF, S ,8, DELTA 

2  FORMAT ( 71H1  2. OEG. EQUAL- RIPPLE  90DEG.DIFF.NW.  HORST  EFFECT  OF 

1COMP.VAR1ATIONS/25HONOMINAL  BANDWIDTH  F2/F 1«»F8.4/25H  NON.  DERIVED 
2  PARAN.  S*»F8.6,5H  3*,F8.6,9F  DELTA*, F4. 21 

3  REA0(5,11  C 

DO  4  1*1,2 

RU  1*I-S^2./SQRT(C(  1 1 l  +  SQRTI S*S/4./C ( I  l+i.-l./C( 1)1 1**2 

4  A(I)*0.5/(1./C( I)+R(I)»-l. 

WRITE(6,5)U( II ,C( I 1,RI II, I* 1,21 

5  FORMAT (25H0N0M.  REALIZ.  CHAN . 1  A*,F3.5,5H  C*,F8.5,5H  R-.F8.5 

1/ 14X,6HCHAN.2,3F13.5/IX, 58(  1H*11 

DO  6  N*l»8 

6  DIM-0. 

GO  TO  12 

7  READI5, 8TI0IN 1 ,N-1 , 71 

8  FORMAT I7F5.21 

IF1DI11.8Q.0. .AND. 0(21. EQ.O. • AND.DI 3 i .EQ.O. • AND.D! 41 *£<3.0.  .AND. 
10(5)*EQ.0.»AN0.D(6)  . EQ.O. « AND.OI  7 1  .Ev'.C.  1  GO  TO  3 
00  10  Ml, 5 
PEPdl-O. 

10  PERI  1 1-PI 

12  DO  13  N-1,7 

I AINl-l 

13  IFIDIN1.NE.0.1  J  AIN1-2 

DO  50  Nl-l, 1 1 

VI 1, 11* l. *011 1/100. *(-l»l**Nl 
VI  1,21*1.-01 l 1/100 »♦(-!. 1**N. I 
DO  50  N2*1,I2 

Y(2»l)*l. +0(21/1 00 »*(-!. )**N2 
VI 2,21*1.-01 2 1/100.*(~1.)**N2 
DO  90  N3*l , 1 3 

Y(3, 11*1. +0(31/1 00 .*<—1. 3**N3 
Y(  3, 21*1. -0(  31/100. *(-l.I**N3 
DO  50  N4*l ,  1 4 

V(4, 11*1. +0(4 1/100. *(-l.l**N4 
Y(4«2)*l.-0(4)/l00.*(-l»  >**N4 
00  5S  N5*i,I5 

Y45tll-l.+0(?l/100.*(-i. )**N5 
Y(5,2l*l.-0!5i/lOO.*(-l.l**N5 
DO  50  NS-1,16 

Y(6»ll*i. +0(61/1 00. *(-i. Sf**N6 
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*TC***5TO*«K»^**«?^^ 


PROGRAM  LISTING  (Continued) 


V<6,2)-l.-D46)/i00.*<-i.)**N6 
00  50  N7-1.I7 

Y47»l>*l.+D<7)/100.*4-1. )**N7 
Y47,2}*l.-0(7i/l00.*t-l« )**H7 
00  15  Ml, 2 

Ell,Il*SQRT4CtI)/Rt  Ill/Yll,t)/Y(5,I)/Y4  7,n-l./SQRT(R4n*C(m 
l/Y<i,ll/V44,I  )-SQRT(R<n*Cm  »/YI2,I)/Y<5,  I> 
EI2»lMl./Yll,n/Y(2,I)/YI4,n/YI5,n 

Et3,n*SQRT<Rm*CI  m*<l./YI2,I)  +  i./Atn/Y(3,n)/(Y(5,n+YC6,I) 
i/A<imi./soRT<R<:)*cm>/Y<  uii/yu.i  MSORTictn/Rt  m/vii.n/ 

2  <Y(5tIM!Yt6,n/Atm 

EU*n*(l./Y(2»I)  +  l./YC3,n/A<n)/Y{  l«l)/Y (4,11/4  VI 5, 1 J*Y<  6,  11/ 
utm 
N*1 

00  30  X>*1«  101 

IF4N.6G.-4.AND.I.LT..  10.1)  GO  TO  30 
X-X 1*4  X2/X1)**(  FLOAT  .I-D/100.) 

P--ARCT  4-X*B*E4 l,l),E42,l)-( X*8 )**2 )+ARCT4 X*B*E43» 1 ),E44»l)-*4 X*B) 
1**2) ♦ ARC T4-X/8*E ( 1 » 2) » E ( 2, 2) -4X/B) **2)-ARCT I  X/B*E  43,2) ,E 44,21- 
24X/B>**2) 

1FI 2-1 )  20*20.21 

PE41MP 

GO  TO  30 

IF4 P-PO )  26,30,22 

IF4N-2)  23,24,25 

PE42MP 

GO  TO  30 

N-3 

GO  TO  30 
PE44I-P 
GO  TO  30 
IFIN-2)  2 7,28,29 
N*2 

GO  TO  30 

PE(3I«P 

GO  TO  30 

N*4 

PO»P 

PE45MP 

00  35  Ml, 5 

IF4PE41  l.GT.PEPt !))  PEP4IMPEII1 
IPIPEin.LT.PENim  PEH4 1  1-PE4 1 ) 

CONTINUE 

iFiotan  5i,5i,60 
00  52  1-1.5 

PE4II»PEU!/PI*180. 

WRITE 4 6,55)  PE 

FORMAT  4 17X.2HFI, 5X, 5HI.MAX, 4X, 3HMIN.4X, 5H2.MAX, 4X,2HF2/ 
115X,3841H«)/13H  OP  NOM./DEG  .5F8.2) 

0481*1. 

GO  TO  7 
00  62  Ml, 5 
PEP4IMPBP4 1) /PI *180, 

PE*(!MPBM(l>/PI*l8Q. 

WRITE46.65)  4 04 N) ,N*1, 7 ) , PEP, PEN 

FORMAT (69HQC0MPQNENT  VARIATION  <*/-PCT>  Rl  R2  R3  Cl  C 
12  C3  0V/23H  (OPPOSITE  IN  CHAN. 1,21, 6X,42< 1M-I/28X,7F6.2 
3/13HOOP  MAX  /OEG  ,5FE.2/13H  OP  MIN  /OEG  .5FR.2I 
GO  TO  7 
JMO 
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APPENDIX  2 

COMPUTER  PROGRAM  FOR  THE  WORST-CASE  ANALYSIS  OF  THE  DOME  FILTER 
CIRCUIT  WITH  OPERATIONAL  AMPLIFIER 

INPUT  DATA  : 


Data 

Format 

Vfl 

F10.G 

c(chan. 1),  c(ehan. 2) 

2F10.6 

ARr  AR2,  ACr  AC2,  AR3,  AR4 

6F5.2 

blank 

» 

c(chan.  1),  c(cban,  2) 

Comment 


Tolerance  in  percent 
(tolerance  set) 

Any  number  of  tolerance  sets 

*\ 


If  analysis  for  different 
circuit  parameters  is 
desired 

Any  number  of  tolerance 
sets 


*As  many  sets  of  datacards  may  be  used  as  desired. 
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PROGRAM  LISTING 


C  DOME  FILTER  (OIFF.AMPL. ) ,  WORST  CASE  * 

DIMENSION  C(2),R(2),IA(6)«D(7),Y(6,2),E(3f2),  l 

1PEP(5),P£M<5),P£(5)  ] 

EQUIVALENCE  < IA<1> ,  I II  ♦  (I  A<  2  )  , 12 ) ,  <  IA<  3)  *  1 3) ,  1 1  AIM  , 1  4)  »  1 

1<IA(5>.I5),(IA<6),I6>  \ 

REAL  K( 2) 

PI*3. 1415926  '  | 

READ45, 1)  F  | 

1  FORMAT C2F 10. 6)  ; 

XI*  SQRT  < 1 . / F ) 

X2*SQRT  <  F  > 

G*<(Xl*X2)»*2-4, 1/U1+X2-2.) 

S*SCIRT<G*S0RT<G*G*2.*G*SQRn  2.*G-4. )))  ; 

B*SQRT ( l.*S*S/2.-G*SQRT  ( t  l.*S*S/2.-G)**2“l.  ) ) 

OELTA*2.*ATAN<<S*S-2.*S*t8-l./B>-(8-i./BJ**2>/tS*S*2.*S*(3-l./B) 

l-<B>l./BI**2n*180./PI 

WR1TEI6,2)F,S,8,DELTA 

2  FORMAT ( 71H1  2.0EG.EQUAL-R  IPPLE  90DEG.DIFF.NW.  WORST  EFFECT  OF 

1COMP.VAR! ATI0NS/25HQN0HINAL  8ANDWIDTH  F2/F 1*, F8.4/25H  NOW.  DERIVED 
2  PARAM.  S*,F8.6,5H  B*,F8.6,9H  DELTA*, F4. 2) 

3  REA0(5,1)  C 

DO  4  1*1,2 

C*SQRT (C(I))*S/2./(l.+C( I)) 

Rm«IG+SQRTIG*G-l./<l.4Cm  M»**2 

4  Km«2./R<I)+2.*Cm  +  l. 

WRITE(6,5HCU),Kin.Rin,I*l»2) 

5  FORMAT  1 25H0N0M.  REAL  12 .  CHAN. I  C*,F8.5»5H  K*,f=8.5,5H  R*,F8.5 

1/14X,  6HCHAN.2, 3F  13.5/1X»58(  IH*)> 

00  6  N*l,7 

6  D(N)*0. 

GO  TO  12 

7  REA 0(5.8) (0(N)*N*1,6I 

8  FORMAT ( 7F5.2) 

IF(D( 1).EQ.0..AND.D(2). EQ.O* . AND.D( 3).EQ.O. .AND.D(4) .EQ.O. .AND. 
1D(5).EQ.0..AND.0(6) .EQ.O.)  GO  TO  3 
DO  10  1*1,5 

PEP( 1 )*0. 

10  PEM( I )*PI 

12  DO  13  N*1 ,6 

I A(N)*1 

13  IF(D(N).NE.O.)  ’<  A  ( N )  *  2 
DO  50  Nl*l, I l 

Y(l,l)*l.«Om/iOO.*(-l.)**Nl 
Y(l»2)*i«-D(l)/100.*(-i. )**N1 
DO  50  N2*l* 12 

Y(2,l)*l.+D(2)/100.*(-1. )**N2 
Y(2,2)*l.-0(2)/100.*(-l.)*»N2 
DO  50  N3*l ,13 

YI3,i)»l.+0(3)/100.*(-l. )**N3 
Y(3,2)»l.-0(3)/100.*(-l.)**N3 
00  50  N4*l , 14 

V<4fl)*l.+0(4)/100.*(“l. )**N4 
Y(4,2)»l.-0(4)/100.*(-l.)**N4 

00  50  N5*l , 15 

Yi5,l)»l.*0(5>/100.*(-1. J**N5 


50 


15 


20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 


35 

50 

51 

52 

55 


6C 

62 

65 


PROGRAM  LISTING  (Continued) 

Hi5^,*U'0<5,/l00**<'L)**N5 

00  50  N6«1,16 

3‘f»iJ*J*+0(6>/IC0.*(-l.>**N6 
^6J2,*l*"0(6|/100,*(-l. )<MN6 
00  15  1*1,2 

^wrli; ! !«!?;  "(??"!$(!(;;<, 1 1 1 

Et3,n*G»H*SQRT(R,n/Ct  I  >  l/Y C  2^T >/V<  3,  J I 

00  30  1*1,101 

IF(N.EQ.4.AN0.1 .IT. 101 )  GO  TO  30 

G*£(2,i )-(X*B)**2 
H*E (2,2 1-(X/B)**2 

IFU-1>  20,20,21 
PEt 1  )=P 
GO  TO  ?0 

IF<P-P01  26,30,22 
IF(N-2)  23,24,25 
P£f 2)*P 
GO  TO  30 
N*3 

GO  TO  30 

PE(4)*P 

GO  TO  30 

1 F l N-2 )  27,28,29 

N-2 

GO  TO  30 
PE ( 3)*P 
GO  TO  30 
N-4 
PO*P 
PE(5I*P 
00  35  1*1,5 

IFtPEtl I.GT.PEPI  1 ) ) 

IFt P€( l ) . LT.PEM4 I)) 

CONTINUE 

1F(0(71»  51,51,60 
00  52  1*1,5 

P£(n*PE(I)/PI*180. 

W#ITE<6,55»  PE 

Dt  7 1*1. 

GO  TO  7 
00  62  1*1,5 

PEP<l>*PEPm/Pl*i80. 

PEPU)*PEH(  I>/PI*iaO. 

WRITE(6,65H0(N),N*1,6)  ,P£P,  PEN 

FORMAT(63HOCOHPON€NT  VARIATION  (fZ-PCT)  R1  R2  Cl 

?13Honp4**»H/ncrP°«TEJN  CHAN.l,2),6X»35(lh-)/28X»6F6.2/ 
213HOOP  MAX  /DEG  ,5Fe,2/13H  DP  MIN  /DEG  ,5 Ffl.2) 

GO  TO  7 
ENO 


P6P(  U*PEtn 
pEMin-PEtn 


C2 
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APPENDIX  3 

SUBPROGRAM  ARCT  REQUIRED  BY  PROGRAMS  IN  APPENDIXES  I  AND  2 


C  QUADRANT-CORRECT  (+/-I80  DEG)  ARCTGtNUM. *D£N0M. )  SUBPROGRAM 
FUNCTION  ARC T ( X I , XR ) 

P 1*3*1415926 
IFtXR)  91 ,94,90 

90  ARCT=*ATANI  XI/XR  ) 

RETURN 

91  !FIXI)  93,92,92 

92  ARCT  =  AT  AN ( X  I /XR )  +PI 
RETURN 

93  ARCT=ATANiXI/XR)-PI 
RETURN 

94  IF(XI)  96,97,95 

95  ARCT=PI/2. 

RETURN 

96  ARCT=-PI/2. 

RETURN 

97  ARC  T'sQ* 

RETURN 

END 


